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ABSTRACT 

The capacity for aerial maneuvering was likely a major influence on the evolution of flying animals. Here 
we evaluate consequences of paravian morphology for aerial performance (Dudley and Yanoviak, 201 1 ; 
Smith, 1952) by quantifying static stability and control effectiveness of physical models (Evangelista et al., 
2014) for numerous taxa sampled from within the lineage leading to birds (Paraves, Xu et al., 2011; 
Gauthier and Padian, 1985). Results of aerodynamic testing are mapped phylogenetically (IVIaddison 
and IVIaddison, 2010; Zhou and Li, 2010; Li et al., 2010; O'Connor et al., 201 1 ; Cracraft et a!., 2004) to 
examine how maneuvering characteristics correlate with tail shortening, fore- and hindwing elaboration, 
and other morphological features. In the evolution of Paraves we observe shifts from static stability to 
inherently unstable aerial planforms; control effectiveness also migrated from tails to the forewings. These 
shifts suggest that some degree of aerodynamic control and and capacity for maneuvering preceded the 
evolution of strong power stroke. The timing of shifts also suggests features normally considered in light 
of development of a power stroke may play important roles in control. 
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INTRODUCTION 

Regardless of how aerial behavior originates, once airborne an organism must control (Smith, 1952) its 
orientation and position in order to safely navigate the vertical environment (e.g., directed aerial descent, 
Dudley and Yanoviak, 2011). Such abilities are present even in taxa with no obvious morphological 
adaptation for flight (Munk, 201 1; Zeng, 2013; Cardona et al., 201 1; Evangelista et al., 2012); at low 
speeds, such as at the start of a fall or jump, inertial mechanisms (Jusufi et al., 2008; Pittman et al., 2013) 
allow for rolling, pitching, and yawing; as speeds increase (or as appendages grow in area), aerodynamic 
mechanisms of control can be employed. Body and appendage configuration and position affect both 
stability, the tendency to resist perturbations, as well as production of torques and forces for maneuvering 
(control effectiveness). In the four-winged Early Cretaceous Microraptor gui, changes in planform, such 
as alternative reconstruction postures or removal of leg and tail feathers, alter stability and the control 
effectiveness of appendages (Evangelista et al., 2014). Furthermore, appendage function can shift entirely 
according to the aerial environment (e.g. asymmetric wing pronation producing yaw at high glide angle 
versus roll at low glide angle) or even completely reverse function (Evangelista et al., 2014). Such results 
are exciting but are based on a single specimen (Xu et al., 2003). Stronger conclusions can be drawn from 
comparative study of several forms within a phylogenetic context. 

One obvious trend in avian evolution is the transition from long tails and feathered legs in early forms 
(Xu et al., 2003; Hu et al., 2009; Longrich, 2006; Christiansen and Bonde, 2004; Zheng et al., 2013; 
O'Connor et al., 2013; Pittman et al., 2013) to later (including extant) forms for which the skeletal tail has 
fused into a short pygostyle and both asymmetric and symmetric flight feathers are absent from the legs. 
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Functional consequences of this shift for aerial maneuvering remain speculative (Smith, 1953; Beebe, 
1915; Thomas, 1997). Similarly, changes in the pectoral girdle have been assumed to enhance a powered 
downstroke (Gauthier and Padian, 1985; Benton, 2005), but may also have influenced maneuvering by 
shifting the center of body mass (Allen et al., 2013) or in enabling the production of wing asymmetries. 
With the exception of (Huynh et al., 2011), previous aerodynamic studies tend to focus on lift and 
drag coefficients and glide angles and specific postures (Chatterjee and Templin, 2007; Alexander et al., 
2010; Koehl et al., 201 1; Dyke et al., 2013), with maneuvering only considered rarely and in single taxa 
(Longrich, 2006; Hall et al., 2012; Evangelista et al., 2014). 

To examine these patterns and to test adaptive hypotheses (Padian, 2001), we can use model tests to 
quantify the effect of shape on static stability and control effectiveness (Evangelista et al., 2014; Koehl 
et al., 201 1; McCay, 2001), using specimens sampled from early paravian (Xu et al., 201 1) and avialan 
(Gauthier and Padian, 1985) evolution (Zhou and Li, 2010; Li et al, 2010; O'Connor et al, 201 1; Cracraft 
et al., 2004). We focus specifically on static stability and control effectiveness; while lift and drag are 
expected to be important in flight evolution, they have been addressed adequately in previous hterature 
(Evangelista et al, 2014; Dyke et al., 2013; Koehl et al., 2011; Alexander et al., 2010). While the capacity 
to generate aerodynamic forces for weight support (Dial, 2003; Dial et al., 2008; Heers et al., 2014) 
was almost certainly developing, we consider here the ability to control aerial behavior (Dudley and 
Yanoviak, 2011). The presence or absence of stability in the various axes and the control effectiveness 
of the appendages should correlate with changes in major morphological features (shortening of the tail, 
enlargement of the forewings) to enhance aerodynamic performance, however, previous work has not 
yet identified what are the patterns. We hypothesize that stability and control are present early in the 
evolution of flight (Dudley and Yanoviak, 201 1); this can be tested by examining the patterns of stability 
and control. Morphologically, for an organism able to perform aerial control with only a weak power 
stroke, we would expect some surface, though not necessarily large, and perhaps with limited force 
generation due to reduced speed, or feather porosity and flexibility (Nudds and Dyke, 2009, 2010; Heers 
et al., 2011)). We would also predict appendages with sufficient flexibility (Kambic et al., 2014) and 
inertia (Jusufi et al., 2008), in combination with aerodynamic forces, to effect maneuvers (Dudley and 
Yanoviak, 201 1). These are all present at the base of the tree of figure 4. 

Alternatively, both stability and control may have been absent early in the evolution of flight, only 
appearing after a strong and bilaterally symmetric power stroke evolved (Dial, 2003; Dial et al., 2008; 
Heers et al., 2014); or very stable platforms with a lack of effective control mechanisms may be observed. 
For these alternatives, morphologically we might expect skeletal features with large surfaces and strong 
muscle attachments including an early carinate sternum, and shoulders with highly restricted ranges of 
motion oriented to provide a fundamental stroke oriented to gravity with limited degrees of freedom (Dial, 
2003; Dial et al., 2008). We may also expect rigidly fused/fixed surfaces located posteriorly to provide 
large stability to offset an inability to generate control. Other then a semilunate carpal assumed to enable 
a flight stroke (Padian, 2001), these are not observed in the taxa tested here (Benton, 2005; Gatesy and 
Baier, 2005; Zhou and Li, 2010; Turner et al., 2012; O'Connor et al., 2012; O'Connor et al., 2013). Our 
hypothesis could be falsified: it is possible for some appendage motions or positions to ineffective or 
unstable, as in the case of certain sprawled postures and leg movements in Microraptor (Evangelista et al, 
2014, though these were also anatomically infeasible). 

Smith (1952), in outlining the potential importance of static stability and control effectiveness, called 
for measurements to test his assertions: "If the conclusions of this paper are accepted the study of 
the remains of primitive flying animals, and even experimental studies on full-scale models of them, 
will acquire a special importance as throwing new light on the functional evolution of nervous systems 
(Smith, 1952)." Subsequent studies of stability have been limited to computational studies with regard to 
aerodynamic function, etc. (Gatesy and Dial, 1996; Taylor and Thomas, 2002; Thomas and Taylor, 2001). 
Computational estimates are useful, but when checked against models (Evans, 2003) or live animals 
(Clark, 2010) there are occasionally unexpected deviations from aerodynamic theory developed for small 
angles of attack and airplane-like morphologies. Therefore, we measured static stability and control by 
measuring the aerodynamic moments exerted on physical models of different specimens in a wind tunnel, 
including at large angles of attack. 
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MATERIALS AND METHODS 

Model construction 

We constructed models (8 cm snout- vent length, figure 1) of four extant birds, seven fossil paravians 
(Xu et al., 2011), encompassing five avialans (Gauthier and Padian, 1985), Microraptor (Xu et al., 2003) 
and Anchiornis (Hu et al., 2009) (Supplemental Material, figure SI) using 3D printing. Fossils were 
selected to sample phytogenies available in 201 1 (when the work was done), although an eighth paravian, 
Zhongornis (Gao et al., 2008), was later dropped due to questions about its phylogenetic position and 
because the specimen was identified to be a juvenile. Additional information regarding the morphology 
of Microraptor, Jeholornis, and Sapeornis has become available since the measurements were done; 
the effect of the new information is addressed in the discussion. To explore parallel evolution and for 
calibration, we also constructed models of three pterosaurs, two bats, and two artificial test objects (sphere 
and weather vane). Construction methods closely followed those of (Koehl et al., 2011; Evangelista, 
2013; Munk, 2011; Zeng, 2013). Solid models were developed in Blender (The Blender Foundation, 
Amsterdam), closely referencing published photographs of fossils and reconstructions from the literature 
(Hu et al., 2009; Longrich, 2006; Hou et al., 1995; Chiappe et al., 1999; Zhou et al., 2002; Zhou and 
Zhang, 2003b; Xu et at., 2003; Zhou and Zhang, 2003a; Zhou and Li, 2010) and casts of Archaeopteryx to 
match long bone, axial skeleton, and body proportions. Modeling was also guided by Starling dissections, 
preserved specimens, and vertebrate anatomy texts (Benton, 2005; Liem et al., 2000). 

Models were printed using a 3D printer (ProJet HD3000, 3D Systems, Rock Hill, SC), then mounted 
on 26-gauge steel armatures with felt or polymer clay filling in gaps between printed parts where flexibility 
was needed for repositioning. Wings were constructed using methods described in (Koehl et al., 201 1; 
Evangelista et al., 2014); wings were traced from the published reconstructions (Hu et al., 2009; Longrich, 
2006; Hou et al., 1995; Chiappe et al., 1999; Zhou et al., 2002; Zhou and Zhang, 2003b; Xu et al., 2003; 
Zhou and Zhang, 2003a; Zhou and Li, 2010), printed on paper and cut. Monofilament stiffening was 
added along feather rachises and the wings were attached to 26-gauge steel limb and tail armatures (scaled 
from the fossils) using surgical tape (3M, St. Paul, MN). This procedure was found in previous work 
to match more laborious and less repeatable application of manually attached real bird feathers (Koehl 
et al., 2011) for fixed-wing tests. The .STL files used to create the models are available for download to 
researchers wishing to replicate our models. 



Body posture and appendage position 

Fossil paravian models were reconstructed with wings spread and legs extended back (Xu et al., 2004; 
Evangelista et al., 2014) (figures 1 and 2). While alternative postures have been considered (specifically 
in Microraptor, Xu et at., 2003; Chatterjee and Templin, 2007; Davis, 2008; Alexander et at., 2010; Hone 
et al., 2010; Koehl et al., 201 1; Habib et al., 2012; Hall et al., 2012; Evangelista et al., 2014; Dyke et al., 
2013; including some now considered infeasible), the aim of this study was to examine maneuvering 
within several species rather than the posture of one, and the legs-back posture is seen in extant birds 
as well as supported for the fossils (Xu et al., 2004; Davis, 2008). For control effectiveness, we tested 
fixed static appendage movements previously identified as being aerodynamically effective (Evangelista 
et al., 2014; Evangelista, 2013): asymmetric wing pronation and supination, wing tucking, symmetric 
wing protraction and retraction, and dorsoventral and lateral movements of the tail (figure 2). The angular 
extent of each movement tested is shown on figure 2. To avoid confusion, we present data for specimens 
as described with all appendages present; artificial manipulations (such as removal of tail and leg surfaces) 
were discussed in (Evangelista et al., 2014; Evangelista, 2013). 

Models were mounted at the estimated center of mass (COM) for the baseline body posture. The 
estimate was formed in Blender assuming a uniform density for the posed model, as in (Allen et al., 
2013). While we did not duplicate the same sensitivity analyses as (Allen et al., 2013), we recognize 
that the COM estimate could vary up to 3-5% of the body length, or by smaller amounts for the variable 
load associated with appendage movements; this uncertainty is usually within the bounds of coefficient 
estimates identified as marginally stable. Mass and wing loading of the organism, while important for 
estimating speed from lift and drag coefficients, do not directly affect the nondimensional coefficients 
presented here. 
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Wind tunnel testing 

Wind tunnel testing used previous methods (Evangelista et al., 2014), with a six-axis sensor (Nanol7, 
ATI, Apex, NC) mounted to a 0.5 inch (12.7 mm) damped sting exiting the model downwind at the 
center of mass (figure 1). In some measurements, a 2 mm steel extension rod or a 3 mm acrylic 
plate were used to avoid geometric interferences and to keep the sting several diameters away and 
downstream of aerodynamic surfaces. The sensor was zeroed for each measurement, eliminating model 
deadweight effects. Models were tested in an open-circuit Eiffel-type wind tunnel with an 18 x 18 x 36-inch 
(45.7x45.7x91.4 cm) working section (Engineering Laboratory Design, Lake City, MN). Testing at 6 
m s^' resulted in a Reynolds number of ^32,000 for all models, matching full scale for Anchiomis, 
Archaeopteryx, Sapeornis, Zhongjianornis, and Confuciusornis. 

Under the test conditions, the aerodynamic coefficients of interest are reasonably constant with 
Reynolds number. Re — UL/v, where L here is the snout-vent length and v is the kinematic viscosity of 
air (Evangelista et al., 2014; Evangelista, 2013). Early in the evolution of animal flight, organisms likely 
flew at moderate speeds and high angles of attack (Evangelista et al., 2014; Dyke et al., 2013) where 
flows appear like bluff body turbulent flows (in which coefficients are largely independent of Re, for 
10^ < Re < 10^). In previous work (Koehl et al., 201 1; Evangelista et al., 2014), we performed a sweep 
of wind tunnel speed, to examine Re from 30,000 to 70,000, to validate that scale effects were not present. 
As additional support for this approach, tests for maneuvering bodies are nearly always tested at well 
below full scale Re, e.g. the largest US Navy freely-maneuvering model tests are well below j -scale. Our 
methods were also previously benchmarked using model tests at full scale Re of gliding frogs (Emerson 
et al., 1990; McCay, 2001) (repeated for comparison), Draco lizards, Anna's Hummingbirds in glide and 
extreme dive pullout maneuvers, hummingbird body shapes in hovering (Sapir and Dudley, 2012), and 
reduced-scale tests of human skydivers compared to actual data (Cardona et al., 201 1; Evangelista et al., 
2012); while at Re ^ 1000, our modeling methods have been benchmarked against extant winged seeds. 
Perching aerial robots, developed to test control algorithms, have shown good agreement between fully 
3D robots and flat plate models with the same planform (Roberts et al., 2009; Hoburg and Tedrake, 2009; 
Tangier and Kucurek, 2005). Results (Evangelista et al., 2014) for lift and drag coefficients using our 
method agreed with those for full-scale Microraptor models in the other modeling tests (Dyke et al., 2013; 
Alexander et al., 2010); our Jeholomis model was at lower Re than Microraptor and is of similar shape. 

Sensor readings were recorded at 1000 Hz using a data acquisition card (National Instruments, 
Austin, TX) (Evangelista et al., 2014). The sting was mounted to a servo (Hitec USA, Poway, CA) 
interfaced to a data acquisition computer, using an Arduino microcontroller (SparkFun, Boulder, CO) and 
specially written code in Python and R (R Core Team, 2014), to automate positioning and measurement of 
windspeed and force/torque. Raw measurements were rotated to a frame aligned with the wind tunnel and 
flow using the combined roll, pitch, and yaw angles by multiplication with three Euler rotation matrices; 
translation from the sensor to the model COM was also included. Transformed measurements were 
averaged over a one-minute recording. We then computed non-dimensional force and moment coefficients, 
static stability coefficients, and control effectiveness (McCay, 2001; Evangelista et al., 2014; McCormick, 
1995). Three series, varying pitch, roll, and yaw, were conducted at 5° increments. Using the automatic 
sting, we obtained 13,792 measurements, with at least five replicates for 18 models in 247 total positions: 
generally 5 each in pitch (88 total), 2 each in roll for two angles of attack (69 total), and 3 each in yaw for 
two angles of attack (92 total). Test positions are indicated in figure 2. 

Static stability was measured by examining the sign of the slope dCm/da (positive slope is unstable, 
negative stable, zero marginally stable, see figure Ic) of the non-dimensional pitching moment coefficient 
C,n near fixed points as the body was subjected to small deflections da (Evangelista et al., 2014; McCay, 
2001; McCormick, 1995): 

pitching moment M = Q.5pU^C,„XS (1) 

where U is tunnel speed, X is the snout-vent length, and S is planform area. Control effectiveness 
(dCm/dd, (Etkin and Reid, 1996; McCay, 2001; Evangelista et al., 2014)) was measured by deflecting 
appendages (figure 2) by an amount d5 and examining the change in pitching moment coefficient. 
Both are unitless (rad^'). Roll and yaw were similarly calculated from series of varying roll angles or 
headings, respectively, with the static stability and control effectiveness partial derivatives taken for the 
roU {0.5pU^C,-XS) and yaw (O.SpU^CyXS) moment coefficients. 

A first-order estimate of maneuvering is obtained from considering the two together and a biome- 
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chanical trade-off is apparent: a stable object can resist perturbations from the environment with minimal 
control effort but will also have difficulty in changing direction (which might be necessary to accomplish 
aerial righting, navigate in cluttered forests, seek resources or avoid predators) (Smith, 1952; Dudley and 
Yanoviak, 2011; Taylor and Thomas, 2002). The metrics underestimate maneuvering in very dynamic 
cases (high advance ratio flapping or where second-order damping terms become important; Sachs, 2005), 
but are adequate for quasi-static maneuvers. Locomotion is a complex task, and passive stability is often 
exploited where possible to reduce control effort; conversely, passive instability may be exploited in 
extreme (and likely elective) maneuvers. The absence of stability, coupled with the presence of large 
control effectiveness, could be used to infer the presence of strong closed-loop neuromuscular control. 
The absence of control effectiveness suggests a lack of control, as even with feedback an ineffective 
surface cannot generate the necessary forces and torques. Thus, while the full control abilities of an 
extinct form are difficult if not impossible to fully enumerate, the simple metrics here provide a useful 
proxy. 

Phylogenetic comparisons 

A Nexus file without branch lengths, available in the Supplemental Material, was assembled from 
pubUshed phytogenies (Zhou and Li, 2010; Li et at., 2010; O'Connor et al., 2011; Cracraft et al., 2004) of 
the study taxa. While revisions to the phylogenetic relationships have been discussed (Xu et al., 201 1; 
Godefroit et al., 2013; Turner et al., 2012; O'Connor et al., 2013), they do not appear to alter the patterns 
in stability and control effectiveness; trees from (Xu et al., 201 1 ; Godefroit et al., 2013; Turner et al., 2012; 
O'Connor et al., 2013) are included in the Supplemental Material. Mapping of discrete maneuvering traits 
as outlined in Padian (2001) was performed in Mesquite (Maddison and Maddison, 2010) with the built-in 
ancestral state reconstruction routines using unordered parsimony. Aerodynamic measurements were 
coded into a matrix giving eight discretized stability values (stable, marginal, unstable, coded according 
to if the 95% confidence interval of measurements includes zero or not) and 12 discretized control 
effectiveness values. The discretized control effectiveness values were obtained from the measurements by 
thresholding based on the moment necessary to overcome measured weather vane stability (dC/d8 > 0.09 
was coded as effective; < 0.09 coded as ineffective), or equivalently, to cause a displacement of the center 
of aerodynamic pressure of about 10% of total length. 

RESULTS 

Representative aerodynamic measurements for pitching stability and control effectiveness are given in 
figure 3 for six paravians and one pterosaur. Tables of all aerodynamic measurements are provided 
in the Supplemental Material. All aerodynamic measurements were coded into both discretized and 
continuous character matrices (see Supplemental Material), which were then mapped onto a phytogeny 
(assembled from Zhou and Li, 2010; Li et at., 2010; O'Connor et at., 2011; Cracraft et at., 2004) to 
examine the evolution of static stability and control effectiveness. The discretized character states for 
pitch, roll, and yaw are shown in figures 4-6. For these figures, color indicates character while the alpha 
transparency (how faded or solid the line appears) indicates the character state. The trees are shown 
to explicitly illustrate the patterns of aerial maneuvering for several different character states within a 
phytogeny (Padian, 2001). Alternate phytogenies (Xu et at., 201 1; Godefroit et at., 2013; Turner et at., 
2012; O'Connor et at., 2013) are given in the Supplemental Material, but do not alter the patterns seen. 
Each trait is subject to uncertainties of measurement (see Supplemental Material), the limitations of 
ancestral state reconstructions with unordered parsimony, and uncertainty in the phytogenies, however, 
in the aggregate the results show consilience (later taxa tines are solid) among pitch and yaw traits as 
discussed further below. 

DISCUSSION 

Additional information and new reconstructions 

The aerodynamic tests described here were performed in the fall of 201 1; since then, more information 
has become available leading to more detailed reconstructions and phytogenies being proposed (Li et at., 
2012; O'Connor et at., 2012; O'Connor et at., 2013; Zheng et at., 2013; Foth et at., 2014). While we are 
not able to perform new measurements, we can theorize about the effect of these newer reconstructions. 
Foth et at. (2014) confirmed Longrich's (2006) finding of teg feathers and does not alter our results. Li 
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et al. (2012) provide a detailed tail reconstruction of Microraptor tail plumage including a graduated 
tail shape with long, midline feathers; we estimate that this tail morphology may have slightly different 
values for stability and control effectiveness but the overall presence/absence pattern we observed here 
would be unchanged. O'Connor et al. (2012); O'Connor et al. (2013) provide further information on 
plumage for Jeholornis; while our forewing reconstruction appears adequate, we lack the proximal fan of 
the "two-tailed" morphology identified in newer specimens. The additional contribution of the proximal 
fan is hard to estimate, and should be considered in future work; however, its forward position suggests a 
small moment arm and therefore a small effect. This is further supported by the marginal stability and 
small control effectiveness of unrealistically sprawled legs, which were also near the center of mass, 
observed in Microraptor models (Evangelista et al., 2014). Zheng et al. (2013) identify distal leg feathers 
and a long/broad tail in new specimens of Sapeornis, while Turner et al. (2012) revise the position of 
Sapeornis to be more basal than Jeholornis. As tested, a more basal position for Sapeornis complicates the 
interpretation of our findings, thus we include a mapping onto such a tree in the Supplemental Material. 
However, taken together, the additional leg and tail plumage described in (Zheng et al., 2013) and the 
more basal position proposed in (Turner et al., 2012; O'Connor et al., 2013) would maintain the patterns 
we see here, shifting one node earlier in the case of pitch (from node 2 to node 1). 

Patterns in longitudinal stability and control 

Long-tailed taxa (figure 3a) show a stable equilibrium point and the tail is effective in generating pitching 
moments, whereas short-tailed taxa (figure 3b) were unstable and had reduced control effectiveness of 
the tail. Notably, the same pattem (i.e., downward sloping C„, versus a) is seen consistently in two early 
representatives of the Avialae, in a long-tailed pterosaur {Rhamphorhynchus), and in the paravian Early 
Cretaceous dromaeosaur Microraptor, suggesting that these patterns more likely derive from shape and 
aerodynamics than from immediate ancestry. 

The study taxa show progressive tail loss as well as loss of leg-associated control surfaces along with 
a concomitant increase in forewing size. Changes in stability and control effectiveness here (as well 
as manipulations in which appendage surfaces were removed with all else held constant (Evangelista 
et al., 2014)) reflect these morphological changes. In pitch (figure 4), taxa shift from being statically 
stable ancestrally to subsequently being unstable (and thus requiring active control, or possibly damping 
from flapping counter-torques; figure 4 red line). Control effectiveness concomitantly migrates from the 
ancestrally large and feathered tail (orange line) and legs to the increasingly capable forewings (yellow 
line), which become relatively larger, gain larger muscle attachments and gain skeletal features and 
stiffness proximally (Benton, 2005) that would improve production of left-right and fore-aft kinematic 
asymmetries needed for control. Distally, bone loss, joint fusion and use of ligaments to reduce degrees 
of freedom (Benton, 2005) would have enabled mechanical feedback and tuned mechanisms as flapping 
developed, enabling neuromuscular control effort to be focused on dealing with increasing overall flight 
instability and active control. 

Transition to forewing control correlates with a significantly enlarged humeral deltopectoral crest 
(Zhou and Li, 2010) and occurs amid progressive acquisition of a fully "avian" shoulder morphology 
(Turner et al., 2012). In addition, the sternum is changing from ossified m Microraptor and Anchiornis 
through varying degrees of loss or ossification without fusion (around nodes 1-2), to ossification, with 
later fusion (node 3) and development of a carinate keel (node 4). Concomitantly, the tail becomes 
much reduced into a pygostyle (Figure 4, node 2) with increased mechanical stiffness (Pittman et al., 
2013), which, combined would have decreased the moments the tail could produce and eliminated inertial 
mechanisms. Other synapomorphies appear at node 4 (figure 4), including a strut-like coracoid and 
triosseal canal (Benton, 2005). Whereas the latter features (node 4) feature in power production, the 
timing of the former features (nodes 1-2) appears consistent with enhanced forewing control effectiveness. 
Ontogenetic tests (Evangelista, 2013) show 4-day post hatching Chukar Partridge are capable of extreme 
maneuvers (rolling and pitching 180°) before strong development of a carinate sternum, suggesting this 
interpretation is correct. 

Roll and yaw control at high angle of attack is present early 

In roll (figure 5), taxa were stable at high angles of attack, but either unstable or marginally stable at low 
angles of attack; large asymmetric wing movements (i.e., wing tucking) were always effective in creating 
substantial rolling moments early, well before development of a power stroke. Also, as animals developed 
the ability to fly at lower angles of attack, active control of roll would have become necessary, perhaps via 
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inertial modes of the tail (only available before node 2) or legs, and later augmented with the forewings 
as they became larger and more capable of left-right asymmetry during full force production (carinate 
sternum, node 4). In all channels, the presence of control effectiveness early is contrary to assertions 
that aerodynamic functions in early paravians are unimportant (Foth et al., 2014). Wing movements with 
high control effectiveness change during evolution in a manner consistent with predicted shoulder joint 
mobility (criterion 1 of Gatesy and Baier, 2005). 

In yaw, most taxa at high angles of attack (figure 6a, green lines) were marginally stable as might 
be expected from symmetry. Taxa with long tails were also stable at low angle of attack (figure 6b, 
purple line), in agreement with computational predictions for similar shapes (Sachs, 2007), but as tails are 
reduced, yaw stability becomes marginal and control migrates from the tail (violet line) to the wings (as 
in pitch control, gray line). Asymmetric wing pronation and supination (figures 6a and 6b) was effective 
in generating yawing moments in all taxa, suggesting maneuverability in yaw early in bird evolution. As 
the tail becomes shorter, flight becomes marginally stable or unstable (node 3), and control effectiveness 
must migrate (as in pitch control) from the shortening tail to the enlarging forewings. In yaw as in roll, it 
is possible that a carinate sternum (node 4) enables more capable left-right asymmetry during full force 
production in extant birds. 

Stability in roll and in yaw shifts with angle of attack, which may serve as a proxy for glide angle or 
the steepness of descent. At high angles of attack, roll is the more stable axis, whereas stability in yaw is 
greater at low angles of attack. Forewing asymmetrical movements created substantial yawing moments 
in all taxa, whereas forewing symmetrical movements were only effective in later taxa. The presence 
of such stability shifts is more consistent with a transition from high glide angles to lower glide angles, 
as predicted by an aerial maneuvering hypothesis (Dudley and Yanoviak, 2011); they are inconsistent 
with a fundamental wing stroke with fixed orientation to gravity. Increased stiffness of the shortening tail 
(Pittman et al., 2013) would also enable high force production at occasional, critical moments of high 
angle of attack flight, such as during landing. The stability and control effectiveness patterns we observed 
illustrate a control system that works at steep angles in ancestral taxa, shifting to one optimized for low 
angles in derived taxa including extant birds. 

Maneuvering and the evolution of flight 

The findings suggest that the capacity for maneuvering characterized the early stages of flight evolution 
(Dudley and Yanoviak, 2011), before forewings with a power stroke fully evolved. Although early 
paravians may have been limited to tight coupling of vertebral and retricial movement (Gatesy and Dial, 
1996), overall gross movement of the large tail of early paravians yielded high aerodynamic control 
effectiveness and the body possessed some degree of stability. Combined with likely dynamic forces 
and torques generated by either tail whipping at the root (Jusufi et al., 2008; Pittman et al., 2013) or 
mild asymmetric or symmetric forewing flapping (flapping limited by less robust skeletal and feather 
morphology or porosity), this suggests that the ancestral organisms were still capable of controlled aerial 
behaviors at high angles of attack (figures 4-6). Gradual evolution of improved maneuvering ability 
(increased control effectiveness, reduced stability) via consistent aerodynamic mechanisms is consistent 
with a continuum of aerial behaviors ranging to full aerial (Dudley and Yanoviak, 20 11; criteria 3 and 4 of 
Gatesy and Baier, 2005). The staggered acquisition of certain morphological characters (e.g. sternum 
ossification; pygostyle) is consistent with aerial maneuvering evolving in an incremental, rather than 
quantum, manner Subsequent shifts in control would be consistent with more shallow glides facilitated by 
incipient wing flapping, which may have served initially in control but then ultimately became the power 
stroke characteristic of modern birds. Incipient flapping may thus have become elaborated as a control 
response (Smith, 1952) to instabilities (figure 4 node 1; figure 6a node 3; figure 6b node 1) demonstrated 
here. Body center of mass was migrating forward (Allen et al., 2013), but this is coupled with loss of 
large posterior surfaces (long tail and leg flight feathers) and coincidence of the wing center with the 
COM. Active control was thus required as static stability was reduced and eventually lost, and associated 
forewing movements would also have enhanced aerodynamic force production and provided a means for 
inertial attitude adjustment. Once the transition to wing-mediated maneuverability and control began, 
larger surfaces and increased musculature would have facilitated dynamic force production for weight 
offset via the power stroke characteristic of modem birds. 
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Figure 1 . Model construction, testing, and measurement of moments. Models were developed in 
Blender (a) from fossils {Archaeopteryx shown) and constructed and tested (b) using previous methods 
(McCay, 2001; Koehl et al., 2011; Evangelista et al., 2014). For simple cases such as a sphere or a 
weather vane, the relationship between slope and stability (c) is observed by plotting pitching moments 
versus angle of attack. Moments for sphere are not statistically different than zero, indicating marginal 
stability as expected, further validating the methods. Models for fossil paravians studied are shown in (d) 
Anchiornis hind limbs rotated out of test position to show plumage for illustration only. 
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Figure 2. Appendage movements tested to determine control effectiveness. Light gray indicates 
baseline posture, dark gray indicates appendage deflection. Appendage movements were selected based 
on those observed to be effective in previous work (Evangelista et al., 2014), including (a) symmetric 
wing protraction (e.g. wing sweep to ±45°); (b) tail dorsiflexion to ±15°; (c) tucking of one wing; (d) tail 
lateral flexion to 30°; and (e) asymmetric wing pronation/supination to 15° (e.g. left wing pitched down, 
right wing pitched up). 
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Figure 3. Representative aerodynamic measurements for pitching stability and control effectiveness. 
Long-tailed taxa (a) have a stable equilibrium point at 10-25° (yellow line) and the tail is effective in 
generating pitching moments at low angles of attack (pale yellow box indicates measurable moments for 
given tail deflections). In short-tailed taxa (b), including extant Lams, the equilibrium point at 0-5° is 
unstable (red line) and the tail control effectiveness is reduced (no measurable moments for the given tail 
deflections). One example (Rhamphorhynchus) drawn from pterosaurs illustrates similar possibilities in a 
phylogenetically distant taxon. 
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Figure 4. Evolution of pitch stability and control effectiveness. Pitching stability is plotted in red hues, 
indicating stable (pale), marginally stable (medium), and unstable (solid). Control effectiveness of the tail 
in generating pitching moments is plotted in orange hues, indicating large control effectiveness (pale) or 
reduced control effectiveness (solid). Control effectiveness of symmetric wing protraction/retraction is 
plotted in yellow hues indicating large (pale) or reduced (solid). Consilience among the three traits (nodes 
1-2) indicates that early in the evolution of Paraves, taxa are stable with a large degree of pitch control 
from the tail; later taxa are unstable, and control has migrated from the now reduced tail, to the wings, 
which become larger and develop skeletal features that would enhance control and the production of 
left-right and fore-aft asymmetries. Nodes 1-4 are discussed further in the text. 




wing tucking: effective, moderate 

Figure 5. Evolution of roll stability and control effectiveness. Characters shown are stability at low 
(15°) angle of attack (mostly unstable due to symmetry; Sapeornis marginal); stability at high (75°) 
angles of attack (all stable); and control effectiveness of asymmetric wing tucking in roll (always 
effective). As animals developed the ability to fly at reduced body angles of attack, more active control of 
roll would have been necessary, at first perhaps using inertial modes of the tail, but subsequently 
augmented with the forewings. 
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Figure 6. Evolution of yaw stability and control effectiveness. At high (75°) angles of attack (a), taxa 
are mostly marginally stable as might be expected for high angles (e.g. at 90° angle of attack all forms are 
marginal). Asymmetric pronation and supination of the wings are always effective in generating yaw at 
high angles of attack. At low (15°) angles of attack (b), by contrast, long-tailed taxa are stable and can 
control yaw with the tail. As tails reduce in size (nodes 1-2), taxa become unstable in yaw at low angles 
of attack and lose the ability to control yaw with the tail. However, asymmetric movements of the wings 
are effective in producing yaw throughout the evolution of this clade, and control would thus have shifted 
from the tail to the forewings, paralleling the shifts seen in pitch. 
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